Cosmo-Particle Searches for Supersymmetry at the Collider Detector at Fermilab 
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Some theories of particle physics are so compelling that it is worth doing a comprehensive and 
systematic set of experimental searches to see if they are realized in nature. Supersymmetry is one 
such theory. This review focuses on the motivation for a broad set of cosmology-inspired search 
strategies at the Tevatron and on their implementation and results at the Collider Detector at 
Fermilab (CDF) with the first few fb _1 of integrated luminosity of data. 



I. INTRODUCTION 

Several theories are candidates for describing particle 
physics which goes beyond the Standard Model (SM). 
Some of them are so compelling that it is worth doing 
a comprehensive and systematic set of searches to see if 
they are realized in nature. Supersymmetry (SUSY) is 
one such theory [1]. In this paper we touch on some of 
the theoretical motivations for SUSY, and concentrate on 
important possible connections to dark matter, astron- 
omy, and cosmology. If these connections are realized in 
nature then cosmology-inspired models, and their impli- 
cations, will help us in deciding where to focus the bulk 
of our experimental search efforts to discover SUSY. This 
paper describes a broad set of cosmo-particle searches for 
SUSY at the Tevatron using the Collider Detector at Fer- 
milab (CDF) [2] . A review covering comparable searches 
at the D0 detector (CDF's sister detector [3]) can be 
found in Ref. [4]. These searches also have important 
implications for when the Large Hadron Collider (LHC) 
at CERN becomes operational. 



A. Supersymmetry and What Makes it 
Theoretically Compelling 

There are many descriptions of SUSY and its moti- 
vations and we refer the reader to excellent reviews [1]. 
We quickly summarize the salient issues for our purposes. 
Essentially, SUSY is a theory that predicts a symmetry 
between fermions and bosons. The Minimal Supersym- 
metric Standard Model (MSSM) describes a scenario in 
which all the quarks and leptons have bosonic counter- 
parts, denoted as squarks and sleptons, q and I respec- 
tively. Likewise, the gauge bosons, including the as-yet 
undiscovered Higgs, have fermionic counterparts, which 
are, as a group, known as the gauginos. The ElectroWeak 
(EWK) eigenstates of the gauginos mix and their mass 
eigenstates are referred to as charginos and neutralinos, 
for the charged and neutral supersymmetric gauge parti- 
cles respectively. We typically order them via their mass 
hierarchy: x?, X°, xt Xa and xf, an d xt ■ 

SUSY helps solve a number of problems in the SM. 
For example, the SM predicts a divergent value for cor- 



rections to the Higgs boson's mass, but Supersymmetry 
offers a way around this problem [5]. When the massive 
sparticles are included in the loops, having bosonic and 
fermionic counterparts make the value of the corrections 
to the mass mostly cancel, leaving behind a finite value if 
the sparticle masses are at or below the TeV scale. An- 
other exciting possibility is that SUSY provides a way 
for the different forces' coupling constants to unify at 
the Grand Unified Theory (GUT) scale, provided that 
the sparticles that have masses below the Fermi scale. 
There is no a priori requirement that this happen, but 
the unification of the EWK and Strong forces at higher 
energies is very suggestive [6] . 

One thing that makes SUSY difficult to look for experi- 
mentally is that it has many different versions with vastly 
different phenomenological features. In particular, as in 
the SM, the particle masses and coupling are not specified 
from first principles. The most model independent search 
method is to pick a set of trial masses and couplings and 
search for each version of the sparticles that way. This 
is not a practical approach experimentally, though, as 
general SUSY models have 100 free parameters [1]. A 
better avenue incorporates theoretical principles and ex- 
perimental constraints, and there are some good general 
phenomenological studies in this vein, for example see [7] . 
One of the most important constraints is that since SM 
particles can decay though other particles in loops, any 
version of SUSY which predicts or allows a short proton 
lifetime (< 10 31 ~ 33 years [8]), must be incorrect. Also, 
any model which predicts a particle's mass to be equal to 
the mass of their corresponding sparticle, as is in the rela- 
tionship between matter and antimatter, must be wrong; 
if sparticles had masses equal to their particle partners, 
we would have discovered bosonic electrons long ago [8]. 
For this reason, it has been argued that SUSY is a broken 
symmetry. 

Identifying the correct methods of symmetry breaking, 
or uniquely predicting the masses and/or couplings has 
proved elusive. One commonly suggested way to protect 
the lifetime of the proton, and solve other problems, is 
the assertion of the conservation of i?-parity [9] , defined 
as R = {— i) 3 ( B -- L )+ 2S ' j where B is the baryon number, L 
is the lepton number, and S is the spin. This gives R = 1 
for SM model particles and R = — 1 for MSSM particles. 
This protects the proton's lifetime because there are no 
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lighter SM particles for it to decay to that do not vio- 
late other symmetries or conservation laws. i?-parity vi- 
olating terms would have to be small for lepton number 
violation and still allow neutrino mixing (lepton flavor 
violation) [10]. A side benefit of i?-parity conservation 
is that it also provides an exciting possibility that SUSY 
could be easily and directly tied to dark matter and cos- 
mology. 



B. The Dark Matter in the Early Universe and 
Today 

1. Astronomy and Dark Matter 

To briefly review the "dark matter problem," we note 
that models of galaxy rotation with stars and hydrogen 
gas have long been known to significantly underpredict 
the rotation velocity of stars at large radii from the cen- 
ter of a galaxy [11]. Simulations that include lots of other 
massive particles, whose dominant interaction is gravita- 
tional, do a much better job of matching the observed 
results much more closely. Since these massive particles 
aren't seen, i.e., the particles don't interact with photons 
or other SM matter much, they are called "dark matter." 
A particle solution would also explain why most of the 
mass is clumped at the center of the galaxy (due to grav- 
ity), but can also form a "halo" around the entire galaxy 
as would be expected if dark matter were a particle that 
has only small coupling which we can't see, such as a 
SUSY particle. 

A second easily understood piece of evidence that the 
dark matter in galaxies is a particle comes from observa- 
tions of collisions of galaxy clusters [12]. Figure 1 shows 
the collision remnants where the mass of the baryonic and 
dark matter components are shown separately from X- 
ray and gravitational lensing methods respectively. The 
two cluster components are well separated after the col- 
lision. The particle physics explanation is that the dark 
matter particles move through each other since they 
rarely interact, while the baryonic matter slows down 
significantly, due to SM interactions. 



2. Dark Matter, Particle Physics, and Cosmology 

Many versions of SUSY predict massive, neutral, sta- 
ble particles that could explain the data. For example, 
if i?-Parity is conserved then, like the proton or the elec- 
tron, the lightest SUSY particle (LSP) cannot decay due 
to conservation of energy. While many models of new 
physics predict dark matter candidates one of the things 
that makes SUSY special is that it provides a well studied 
model of particle physics that can also make quantitative 
predictions for astronomy and cosmology. This includes 
predictions about early universe physics, including the 
baryon asymmetry as well as a prediction for the dark 




FIG. 1: An image of the dark matter and baryonic mass dis- 
tributions after the collision between two clusters of galaxies. 
The blue represents the amount of matter present, as mea- 
sured by gravitational lensing, but not otherwise seen. This 
has been attributed to dark matter. The red represents the 
mass distribution, as measured by X-rays from the galaxy 
remnants, mostly hydrogen gas. The interpretation of this 
data is that the dark matter has largely passed through un- 
affected by the collision, but the baryonic matter has signifi- 
cantly slowed due to SM interactions. Taken from [13]. 



matter relic density imprinted in the cosmic microwave 
background radiation. 

Cosmology measurements have entered a precision era. 
WMAP and other experiments currently have SM mat- 
ter accounting for about 4% of the mass of the known 
universe, and dark matter accounting for about 23%; the 
remaining 73% of the universe is called dark energy [14] . 
While an analysis of dark energy is outside the scope of 
this review, the dark matter relic density, and the ratio 
of the mass scale of the dark matter and its coupling with 
ordinary matter, provides strong constraints on SUSY if 
indeed the lightest SUSY particle is the dark matter. If 
the dark matter is made of sparticles, there is the excit- 
ing possibility that SUSY not only provides a dark mat- 
ter candidate, but also describes early universe physics 
and can provide full calculation of the dark matter relic 
density, ^susy dm [15]. 



C. Overview: A Number of Different Types of 
Cosmo-Particle Solutions 

If nature has chosen the simple solution that 23% of 
the mass of the universe is a single type of particle that 
is, in fact, the lightest SUSY particle, then this places 
very severe constraints on the types of SUSY that can 
produce this scenario. That being said, there are still 
many different types of solutions possible, even with this 
assumption. We next consider some possibilities, sum- 
marized in Figure 2, that describe some of the crucial 
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FIG. 2: Different decays of SUSY particles to the LSP in R- 
parity conserving scenarios with different lifetimes (r). (a) 
shows the fast decay of a selectron to the lightest neutralino, 
which is a cold dark matter candidate in mSUGRA models, 
(b) shows the decay of a long-lived stop, a CHAMP in this 
case, to the lightest neutralino. (c) shows the decay of the 
lightest neutralino to a photon and a light graviton, which is 
a warm dark matter candidate in GMSB models. Note the 
wide variation of time scales through which these processes 
occur. This has important implications for both the early 
universe physics and the dark matter content of the universe 
today. 



elements of how the early universe might have evolved 
into the universe we know today, and how sparticles, and 
their properties, fit in. 



1. Minimal Solution with Cold Dark Matter 

A minimal scenario is that the lightest sparticle is neu- 
tral, weakly interacting, stable, and mostly produced in 
the early universe (t = 10 -24 sec) as the decay product 
of all the heavier, short-lived sparticles (see Figure 2a). 
It is still here today since it cannot decay and has such 
small interaction cross sections that it is rarely annihi- 
lated. The model is minimal in that all the dynamics 
are mostly over after inflation [16]. Cold dark matter 
is favored by most cosmological models [17], and many 
SUSY models have a lightest particle that fits this de- 
scription in the 100 GeV/c 2 mass range, typically the 
lightest neutralino, Xi [!]■ The minimal model of SUSY 
that also includes potential grand unification with grav- 
ity is known as Minimal Supergravity, or mSUGRA for 



short [18]. This model has become popular in the field 
for these reasons, and because it is well described by only 
a few parameters, described in more detail in Section III, 
we refer to it as our baseline, minimal search model. 



2. Non-Minimal Solution with Cold Dark Matter 

While the minimal models have the advantage of ele- 
gance, there is no compelling reason to believe that na- 
ture has chosen such a simple solution. For example, 
it would imply that the amount of dark matter, which 
is much more prominent in the universe than SM parti- 
cles, has simply evolved along with the universe since 
10~ 9 seconds after the Big Bang, when the tempera- 
ture dropped below levels where the other sparticles were 
typically produced. Normal matter, especially quarks 
and atoms, have very different roles in the universe over 
time. There are many non-minimal models that could 
produce the dark matter we observe today. For exam- 
ple, there are models with long-lived CHArged Massive 
ParticleS (CHAMPS) that decay into dark matter (see 
Figure 2b) [19]. In this case the dominant SUSY parti- 
cle of the universe on the timescale of universal inflation 
is a CHAMP, but it decays into the LSP/dark matter 
after a nanosecond, a second, or even a month. The ad- 
vantage to these models are that they essentially give an 
extra "knob" which helps in fitting models to the data. 
These models have recently gained favor in parts of the 
theoretical community for other reasons as well [19]. 



3. Non-Minimal Solution with Warm Dark Matter 

The lifetime of the LSP is clearly an important param- 
eter of the theory, but the mass of the LSP may be even 
more important. Warm dark matter, a light particle, is 
also consistent with astronomical data and inflation mod- 
els [17]. An example model that has this feature is Gauge 
Mediated Supersymmetry Breaking (GMSB), where the 
process Xi^ lG is allowed [20], where the G is the grav- 
itino, the SUSY partner of the as-yet unobserved me- 
diator of gravity, the graviton. As in the non-minimal 
cold dark matter solution, the Xi is stable on the time 
scale of the early universe but decays on the nanosecond 
timescale to the keV mass G, which is still here today, as 
shown in Figure 2c. 



Other Possibilities, Including the one that Dark Matter 
and SUSY are Unrelated 



Perhaps the situation is far more complicated than we 
thought and dark matter and SUSY are unrelated, but 
SUSY is still a correct description of nature. Maybe there 
are two (or more) massive, stable, weakly interacting par- 
ticles in nature [21]. Possibly Axions play a role [22]. In 
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this case, we must look for the most general models in- 
cluding i?-parity violating terms [9, 10]. These scenarios 
are also outside the scope of this paper. 



D. Outline 

With this context we next discuss a number of cosmo- 
particle searches for sparticles at the Tevatron, as sum- 
marized in Table I. After a quick experimental intro- 
duction in Section II we focus on mSUGRA searches for 
convenience in Section III. These include separate direct 
searches for light squarks and gluinos, gaugino pair pro- 
duction, as well as for Sbottoms and Stops. There are 
also important indirect searches and direct searches for 
CHAMPS, which are are particularly sensitive to certain 
regions of parameter space. In Section IV, we discuss 
GMSB models. 
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Gauge Mediated 


Short-lived x? 


IV A 


[30] 


Symmetry Breaking 


in "fffir 






with x? — ► 7G 


Long-lived x? 
Delayed Photons 


IV B 


[31] 



TABLE I: Topics covered in this review. 



II. EXPERIMENTAL ISSUES 

Before describing the searches we discus some other ex- 
perimental issues in looking for the particles of the early 
universe. The Fermilab Tevatron collides protons and 
anti-protons with y/s = 1.96 TeV and is the high energy 
frontier until the LHC is turned on. Since a tempera- 
ture of 100 GeV corresponds to about 10 ps after the 
Big Bang, the CDF and D0 detectors allow us to look 
back at conditions as far back as a picosecond. As this 
manuscript was created during the Summer of 2009, the 
Tevatron had delivered about 6 fb _1 , the detectors had 



acquired about 5 fb _1 , and analyzed 2-3 fb _1 . This pa- 
per will cover the results from CDF, and point the reader 
to the D0 results for completeness. The Tevatron has re- 
sumed taking data as of September 2009, and the plan 
is to run for at least a year. There are rumors of having 
it run until 2012, as it keeps delivering more and better 
data as time goes on. This, of course, critically depends 
on the LHC's progress; for more on LHC start-up issue, 
see Appendix A. 

A detailed description of the CDF detector can be 
found elsewhere [2]. It is the prototypical collider detec- 
tor with concentric sub-detectors used for particle identi- 
fication and 4-momentum measurement. From the inside 
out it contains a silicon tracker, a central tracker, electro- 
magnetic calorimeters, hadronic calorimeters, and muon 
chambers. After more than 20 years of experience there 
now exists well understood and standardized identifica- 
tion methods for electrons, muons, taus, photons, jets, 
6-jets, and missing energy. The missing energy is par- 
ticularly important for the identification of dark matter 
candidates. As a weakly interacting neutral particle, it 
would, if produced, leave the detector without interact- 
ing, like the SM neutrino, and create the missing energy 
imbalance we call fir- Other new detector techniques 
have come online for Run II, including 25 ps charged par- 
ticle timing [32] , 500 ps photon timing [33] , and powerful 
new fir resolution models [34]. 



III. MINIMAL SUPERGRAVITY MODELS: 
MSUGRA 

Models of SUSY where the symmetry breaking is me- 
diated by the gravity sector are known as Supergravity 
models. Most have some simplifying assumptions and the 
most commonly used minimal model is known as minimal 
Supergravity or mSUGRA [18]. In this model the scalars 
are assumed to have a common mass, mo, and sfermions 
have a common mass, m 1 / 2 , at the unification scale. To 
first order most sparticle masses scale with these parame- 
ters. At incur there are three other free parameters that 
determine the sparticle masses at the EWK scale. They 
are tan/3, which is the ratio of the Higgs Vacuum Expec- 
tation Values (VEV), which will be important later, Aq, 
the common trilinear coupling factor, and sgn(/i), where 
(i is the Higgs mass parameter. While this model may 
be too simple to be true, it provides a useful phenomeno- 
logical benchmark for use in both focusing searches and 
interpreting the sensitivity of the searches in the case of 
null results. 

In a typical mSUGRA scenario available at Teva- 
tron energies, and taking into account experimental con- 
straints, the squarks and gluinos are the heaviest sparti- 
cles, the first and second generation sparticles are mass 
degenerate, and the lightest neutralino is the LSP, mak- 
ing it a good cold dark matter candidate. The typical 
mass relations are approximately m^o : m^o : fn~± : rrig 
= 1 : 2 : 2 : 6 and m,q = nig. The masses of the stop, 
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sbottom, and stau are very dependent on tan/3 and for 
large values, tan/3 > 10, can be very light. This case has 
the aesthetically pleasing aspect that, if true, as in the 
SM, members of the third generation would have very 
different masses than members of the other generations. 

While it's important not to place too much emphasis 
on a particular SUSY model, likelihood fits of mSUGRA 
to existing data including the current Higgs mass limits, 
g — 2 and other experimental data in the plane of m-1/2 
and tan/3 show better fits with high tan/3 [35], as shown 
in Figure 3. This strongly suggests that to have better 
coverage for a variety of models that we do a range of 
searches to cover both low tan/3 and high tan/3. 

Given these inputs our benchmark cold dark matter 
searches at the Tevatron can be broken down into a num- 
ber of types which have distinct advantages and disad- 
vantages, depending on what is realized in nature. They 
are: 

• Direct production of squarks and gluinos: Like 
quarks and gluons, as strongly interacting particles, 
squarks and gluinos should have the largest produc- 
tion cross sections per unit mass at the Tevatron. 
However, as the production cross sections also drop 
rapidly as a function of sparticle mass, they may 
be too heavy to be produced at the Tevatron. An 
additional complication is that since they typically 
produce multijet + fir final states they are difficult 
to separate from the large QCD backgrounds with 
fake fir. 

• Direct production of gauginos: The gauginos are 
typically lighter than the colored objects which 
could make their production cross section large. 
However, they have EWK production mechanisms, 
which make the cross sections smaller per GeV/c 2 



of mass. If the lightest neutralino is about, or 
above, 75 GeV/c 2 (from LEP [36]) then, from 
the above mass relations, the squarks and gluinos 
would be too heavy to be produced at the Teva- 
tron. However, gaugino pairs could be produced at 
observable levels. A second advantage of gaugino 
pair production is that they often have leptonic fi- 
nal states, for which there are smaller backgrounds. 

• Third generation searches: At high tan/3 the stop 
and sbottom masses might be light enough to make 
them easier to produce. Often they produce final 
state b or t quarks which could be identified. Simi- 
larly, gaugino branching fractions to r's can rise to 
100% as the stau gets light making their signature 
also very distinct. On the flip-side, the purity and 
efficiency for r and 6-quark identification are not at 
the same level as their first and second generation 
counterparts. 

• Indirect searches via sparticles in loops: These 
searches measure the branching ratios of SM par- 
ticle decays for evidence of sparticles in loops. If 
all the sparticles are too heavy to be directly pro- 
duced, this may be our best, or only, bet. 

• CHAMPS searches: If the next-to-lightest SUSY 
particle (NLSP) is long-lived, then, as in indirect 
searches, there will be no direct fir signature. In 
this case, this again may be our best, or only, bet. 

We consider all these types separately to maximize cov- 
erage and start with low tan/3, moving to searches with 
higher tan/3. We again note that, as for many CDF re- 
sults, there are typically comparable D0 search results. 
We will not cover these individually, but reference them 
for completeness. 
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FIG. 3: A likelihood fit of mSUGRA in the m 1/2 vs tan/3 
plane to the experimental data. The data clearly favor high 
tan/3 in these types of models. Taken from Ref. [35] 



A. Low Tan/3 Searches 

1. Squarks & Gluinos 

For squark and gluino masses available at the Teva- 
tron, there are three main production diagrams: qq, gg, 
and qg. As the two masses are expected to be similar, 
small differences between them can produce final states 
which can vary significantly. Specifically, the decays are 
expected to be q — > qxl and g — > qq — > q(qXi) where the 
intermediate squark is real or virtual, depending on the 
masses. These produce three separate final states: two 
jets (for nig 3> mq), three jets (for nig = rriq), and four 
jets (for m-g <C niq), where the weakly interacting Xi' s 
leave the detector creating fir- This suggests a multijet 
+ firr analysis. CDF now uses a new, optimized strategy 
with 2.0 fb _1 of data and searches for two jets + fir, three 
jets + fir, and four jets + fir events separately, but with 
a unified analysis to give the best coverage [23] . The anal- 
ysis selects multijet events with large fir and Hr, (where 
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Ht is defined as the scalar sum of the jets and fir) as well 
as restrictive kinematics to ensure well measured fir in 
events. Despite huge QCD production cross sections, the 
sophistication of the analysis techniques have reduced the 
QCD backgrounds with fake fir to the 25% level. The 
rest is ti and other EWK processes with final state neu- 
trinos which produce real jets and fiq<. Figure 4 shows 
the Ht and fir distributions in large Hr events for two, 
three, and four jet events. Table II gives the results of 
the counting experiments. 





Jets 




2 


3 


4 


Selections 


H T > 330 GeV 
fir > 180 GeV 


H T > 330 GeV 
fir > 120 GeV 


H T > 280 GeV 
fir > 90 GeV 


Expected 
SM Events 


16 ± 5 


37 ± 12 


48 ± 17 


Observed 
Events 


18 


38 


45 



case of nig = rriq. If more model dependence is added by 
including the production cross section (i.e., more physics 
in the interaction vertices) mass limits can be extracted 
from where the cross section limits are below the produc- 
tion cross section. This is essentially where the two lines 
in Figure 5 cross, but also takes into account systematic 
uncertainties on the production cross section. The exclu- 
sion region in the squark vs. gluino mass plane is shown 
in Figure 6. Some notes are in order: nig < 280 GeV/c 2 
is always excluded and m > 392 GeV/c 2 is excluded 
when nig — rriq. A typical last step in the process is 
to invoke the mechanism of SUSY breaking to realize the 
sensitivity to model parameters. For mSUGRA, for the 
parameter choice used, the exclusion region is shown in 
Figure 7 in the mo vs. m 1 / 2 plane. Of particular inter- 
est is that the limits go beyond those from LEP [36] for 
mSUGRA models in the region 75 < m < 250 GeV/c 2 
and 130 < m 1 / 2 < 170 GeV/c 2 . Comparable results from 
D0 can be found in Ref. [37]. 



TABLE II: The number of events in the multijet + fir search 
for squarks and gluinos. For example signals see Figure 4. 
Taken from [23]. 

Since there is no evidence for new physics, 95% confi- 
dence level limits are set. Setting limits can be done in 
a number of different ways to illustrate the experimen- 
tal sensitivity to new physics as more and more model 
dependence is added into the results. Cross section lim- 
its have a significant advantage in that they are fairly 
model independent and essentially reflect a balance be- 
tween an acceptance model that comes from the kine- 
matics of the particles produced and their decays, and the 
backgrounds and the ability to reject them. The observed 
and expected cross section upper limits for the squark 
and gluino searches are shown in Figure 5 for the special 
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mSUGRA signal: 
M g = 350 GeV/c 2 
m! = 385 GcV/l- 2 
A ( ?= O, tanp = 5, it<Q 
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H | [GeV] E T [GeV] 

FIG. 4: The H T and fir distributions in (Top) 2 jet + fir, 
(Middle) 3 jet + fir, and (Bottom) 4 jet + fir events along 
with the background expectations and an example distribu- 
tion for squark and gluino production. Taken from [23]. 
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FIG. 5: The 95% confidence level cross section upper limits 
for the case m-g = rriq from the unified multijet + fir search. 
Taken from [23]. 



2. Gaugino Pair Production 

Direct gaugino production at the Tevatron has EWK 
production mechanisms, but can dominate the sparti- 
cle production cross section if the squarks and gluinos 
are too heavy to be produced at the Tevatron. It is 
favored in mSUGRA scenarios where the Higgs mass 
limits and lightest neutralino mass limits from LEP are 
taken into account [38]. The dominant production dia- 
grams are typically XiXi an d X1X2 production, which 
are analogous to SM WW and WZ production. A com- 
mon leptonic decay mode of the chargino is typically 
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FIG. 6: The 95% confidence level exclusion region in the 
squark vs. gluino mass plane from the multijet + fir search. 
Taken from [23]. 
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Most previous Tevatron trilepton + frp searches at 
low tan/3 considered only the eee, eep,, efifi, or fifj.fi 
final states. Rather than focus on e's and fi's sepa- 
rately, the latest search, with 2.0 fb -1 of data, consid- 
ers five separate event quality types in a unified analysis 
to expand the coverage and include final state hadronic 
t's [24]. "Tight" lepton identification helps select elec- 
trons or muons in a manner that produces very pure sam- 
ples (i.e. low backgrounds from jets that fake leptons) at 
the cost of efficiency. "Loose" lepton identification ac- 
cepts more electrons and muons with better efficiency at 
the cost of purity, i.e. has more background from fake 
leptons. A third lepton category is to only require an 
isolated, charged particle track to identify an electron, 
muon or single-prong hadronic tau decay. Note that lep- 
tons considered tight cannot simultaneously be identi- 
fied as loose or a track. Similarly, a lepton identified as 
loose-but-not-tight cannot be identified as a track. The 
five non-overlapping categories, along with the results of 
the background expectations and the experimental re- 
sults are shown in Table III. Figure 8 shows the ipr dis- 
tribution of the 2 Tight + 1 track final state. There is 
no evidence for new physics. 

The interpretation of the data is not only gaugino 
mass dependent, but is also determined by the X2 decay 
branching ratios which are very slepton mass dependent. 
Figure 9 shows the 95% confidence level upper limits on 
the cross section as a function of the chargino mass in the 
large and small slepton mass scenarios. Taking these ef- 
fects into account we find the mSUGRA exclusion region 
shown in Figure 10. The middle region, where there is 
no sensitivity, is where the stau mass changes from being 
above the X2 mass to below. In this region there is such a 
small mass difference between the X2 an d slepton masses 
that there is very little energy for one of the lepton Pt 
in X2 ~* ^ decays, making it undetectable in the detec- 
tor. Reference [39] gives a method for model-independent 
interpretation along with numerical formulae for general- 
izing the Tevatron trilepton search results. Comparable 
D0 search results can be found in [40]. 



FIG. 7: The 95% confidence level exclusion region in the mo 
vs. mi/2 mSUGRA mass parameter plane from the multijet 
+ fr search for squarks and gluinos. Taken from [23]. 



xf — > VF*^? — > (lv)xi, which means that chargino pairs 
will produce a two lepton + frp final state, which is dif- 
ficult to distinguish from SM backgrounds. The domi- 
nant decays of the neutralino are X2 ~ > H ~^ {^Xi) I an d 
X2 —> Z*x1 — » (l + l~)x.i> depending on the m^o — mj 
mass differences. Thus, XiX^i production can result in a 
final state of three leptons + Jpr, which has very few SM 
backgrounds. 



Channel 


Background 


Obs 


3 Tight 


0.49 ± 0.04 ± 0.08 


1 


2 Tight + 1 Loose 


0.25 ± 0.03 ± 0.03 





1 Tight + 2 Loose 


0.14 ± 0.02 ± 0.02 





2 Tight + 1 Track 


3.22 ± 0.48 ± 0.53 


4 


1 Tight + 1 Loose + 1 Track 


2.28 ± 0.47 ± 0.42 


2 



TABLE III: The five final state configurations in the united 
search for chargino-neutralino production and decay in the 
trilepton + Jpr final state. Taken from [24]. 
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FIG. 8: The fir distribution in 2 Tight + 1 Track events as 
part of the search for gaugino pair production in the trilepton 
+ fir final state. Taken from [24] . 
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m„ (GeV/c 2 ) 

FIG. 10: The mSUGRA exclusion region from the search for 
chargino-neutralino production and decay in the trilepton + 
fir final state [24]. Note the unexcluded region in the middle 
which occurs when the stau mass is approximately equal to 
the %2 mass, producing a third lepton which has a Pr which 
is so low as to be undetectable. For more detail, see [39]. 
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FIG. 9: Cross section limits on chargino-neutralino produc- 
tion and decay in the trilepton + fir search for two different 
values of mo. Note that in the bottom plot the limits get less 
restrictive near the m? 



m^o mass point. Taken from [24]. 



mated from data using sideband techniques. Of particu- 
lar note is that the backgrounds have not successfully 
been estimated from Monte Carlo simulations, which 
makes predictions for sensitivity at the LHC precarious. 
An a priori analysis considers three neural net bins, 
shown in Figure 11, and 5 separate mass bins; each of 
which is compared to expectations. The neural net bins 
are defined such that values close to 1.0 are B s — > 
like and values close to 0.0 are background-like. 

Since there is no evidence for new physics a 95% confi- 
dence level upper limit of 5.8 x 10~ 8 is set on the branch- 
ing fraction. Comparable results from D0 can be found 
in [44]. In mSUGRA, cosmology constraints point to high 
tan/3 and the co-annihilation region, where t~x\ mass dif- 
ference is expected to be less than 20 GeV/c 2 /c 2 [41], as 
shown in Figure 12. These limits correspond to an exclu- 
sion of mi/2 below 380 GeV/c 2 and gluino masses below 
925 GeV/c 2 . As expected these are well above the direct 
search mass limits and are the world's most sensitive. 



2. Sbottom Searches 



B. High Tan/3 Searches 

1. Indirect Searches 

The search for B s — > /x + fj,~ is perhaps the most sen- 
sitive to SUSY at high tan/3 since sparticles can show 
up in loops as the branching fraction rises as tan 6 /3 [41]. 
In the SM, the decay of B s — > is heavily sup- 

pressed and has a prediction of BRsm{B s — > M + aO = 
(3.5 ± 0.9) x 10~ 9 [42]. This was the first SUSY search 
out of CDF in Run II [43] and since then has been heav- 
ily optimized using neural net techniques and 2.0 fb _1 of 
data [25]. The backgrounds are combinatorial and esti- 



For high tan/3 the bottom squarks can be lighter than 
the other squarks and may well be the best direct search 
method. Also, powerful 6-tagging techniques, developed 
for top quark properties measurements [45] , enable com- 
plex searches for events from b — > b\i decays producing 
6-jets + fir final state. There are two main types of sbot- 
tom searches. The first is from gluino pair production via 
gg — > (bb)(bb) — > the second is direct sbot- 

tom pair production bb — » (&Xi)(^Xi)- Direct sbottom 
searches in Run II were published in 2007 by CDF [46], 
with an update in 2009 [47] and in 2006 by D0 [48]. 

New results have recently appeared in gluino-mediated 
scenario searches in the multi-6-tag + fir final state [26] 
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FIG. 12: The mSUGRA region for large tan/3 and the exclu- 
sions from various other experiments. Taken from [41]. The 
current limits from Ref. [25] are 5.8 x 10 -8 . 



with 2.5 fb _1 of data. The gluino production searches 
are complementary to the direct sbottom searches which 
are gluino mass independent, but are only sensitive to 
lower mass sbottoms. Two separate optimizations are 
performed to cover the cases when the gluino-sbottom 
mass difference is large, producing larger Et 6-jets in the 
final state, and when the mass differences are smaller and 
the 6-jets are less energetic. After all requirements, in- 
cluding 6-tagging, large Ht, and large fir and neural net 
techniques, the backgrounds are roughly half QCD and 
half top quark and EWK production. The results of the 
two counting experiments are given in Table IV. Since 
there is no evidence for new physics, cross section limits 
on the gluino are shown in Figure 13. Figure 14 shows 
the exclusions in the sbottom mass vs. gluino mass plane 
along with results from Ref. [46] . The most recent results 
in this final state from D0 can be found in Ref. [49]. 



Background 


Large Sm 


Small Sm 


Source 


Optimization 


Optimization 


Electroweak boson 


0.17 ± 0.05 


0.5 ± 0.3 


backgrounds 






Top- quark 


1.9 ± 1.0 


0.6 ± 0.4 


Light-flavor jets 


1.0 ± 0.3 


0.6 ± 0.1 


Heavy flavor 


1.6 ± 0.8 


0.7 ± 0.3 


multijets 






Total expected SM 


4.7 ± 1.5 


2.4 ± 0.8 


Observed 


5 


2 



TABLE IV: The results of the counting experiments in the b- 
jets + tfh; search for gluino-mediated production of sbottoms 
pairs. Large 8m refers to scenarios where mj = 335 GeV/c 2 
and mj = 260 GeV/c 2 , which produces high Et final state 
b-jets. Small Sm refers to scenarios where mgr = 315 GeV/c 2 
and mj = 335 GeV/c 2 , which will produce low Et final state 
jets. Taken from [27] . 



3. Stop Searches 

The lightest stop mass eigenstate, ti, can be sig- 
nificantly lighter than all the other non-LSP sparticles 
due to a large top- Yukawa coupling, especially at large 
tan/3. There have been many searches for pair produc- 
tion and decay of t\ in different decay modes including via 
h -> xtb -> (Zi>Xi)&, h -> x£b -> (jjXi)b, and ex?. 
Searches where both stops decaying t — » c\i were pub- 
lished in [46] , and updated in 2009 [50] . A similar search 
from D0 can be found in Ref. [51]. A search via the 
chargino channels where one stop decays leptonically 
and the other hadronically was carried out by D0 in 
2006 [52]. Recently a new search in the dilepton + jets 
+ ]/<}? final state signature was performed. Special mo- 
tivation for this channel comes from the Run I obser- 
vation that some of the dilepton events didn't "look" 
like top quark pair production and decay [53]. Simi- 



10 



.Q 
Q. 



■a 1 

o 

CD 

CO 

co 
co 
o 

6 10" 1 r 



CDF Run II (2.5 fb" 1 ) 

■ 95% CL limit (m[b]=250 GeV/c 2 ) 
' Expected limit 

95% CL limit (m[b]=300 GeV/c 2 ) 
Expected limit 

g->bb(100%BR) 
b->bx°(100%BR) 




pp^ggat\[s=1.96 GeV 

10" 2 F PROSPINO NLO (CTEQ6M) 

^ R = n F = m(g) 

m(q>500 GeV/c 2 m(x°)=60 GeV/c 2 



_±_ 



250 



300 



350 



400 



m(g) [GeV/c 2 ] 



~ 350 

> 

CD 

^300 

co 

co 

CO 



E 
o 

-I— • 
-t— • 

o 
CO 



250 



200 



150 



100 



7 CDF Run II (2.5 fb 


) y 


g^bb(100%BR) 


^395% CL limit 
Expected limit . 




b^bx°(100%BR) 


: 




■ m(x°) = 60 GeV/c 2 
i m(q) = 500 GeV/c 2 


- 






■/ CDF Run I1 156 pb7 
Excluded Region I 


DO Run II 310 pb" 1 
Sbottom Pair Production 
Excluded Region 


! CDF Run I excluded 

" , i , , , , i , , , , i , , , , i , , , , i , , , 



200 



250 300 350 400 



Gluino Mass [GeV/c 2 ] 



FIG. 13: The 95% confidence level cross section limits for pair 



production of gluinos that decay via g — 
b-jets +tf)r final states. Taken from [27] 
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FIG. 14: The exclusion regions in the sbottom mass vs. gluino 
plane from searches in the b-jets + fhr final state for gluino 
mediated sbottom production. Also shown are results from 
searches for direct sbottom production. Taken from [27]. 



larly, an admixture of stop events among the top quark 
events would effect the kinematics of the final state ob- 
jects and would explain differences in the top mass mea- 
surements in the jets+^r, the lepton+jets+^r, and the 
dilepton+jets+^T final states as seen in early Run II [54]. 

A powerful new method [28] does a fit of the kine- 
matics of all dilepton + dijet + fir events after minimal 
top quark event rejection requirements to search for stop 
events mixed in with the top quark events. In this search 
in 2.7 fb -1 of data the t\, Xii an d Xi masses are allowed 
to vary in the fit. Figure 15 shows the best fit t\ mass 
distribution along with the background predictions. As 
there is no evidence of new physics, branching ratio lim- 
its are set as a function of the t\ and xf masses. The 
results are shown in Figure 16 and are the first exclusions 
on these sparticles in the t\ — > xtb ~^ (^Xi)^ channel. 



C. CHAMPS 

There is a new emphasis in the theory community 
about the role of long-lived sparticles in the early uni- 
verse that decay into the dark matter particles we ob- 
serve today [19]. These particles can be identified in the 
detector because they are heavy, charged, and weakly in- 
teracting. These properties make them interact in the 
detector in ways that make them look like muons that 
are just traveling significantly slower than the speed of 
light. To search for these particles, timing techniques are 
used to measure the time of arrival for muon-like parti- 
cles that arrive later than "expected" [29]. The results 
from 1.0 fb -1 , converted to the "measured mass" are 
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FIG. 15: The mass distribution of dilepton + fr + jet events 
reconstructed under the stop hypothesis. Taken from [28] . 



shown in Figure 17, along with SM expectations. There 
is no evidence for new physics. The 95% confidence level 
cross section upper limits where the CHAMP is a stop are 
shown in Figure 18; similar results are expected when the 
limits are reinterpreted in stau and chargino scenarios. 
Comparable results from D0 can be found in Ref. [55]. 
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FIG. 16: The 95% confidence limits on stop production and 
decay as a function of chargino and ti masses from the search 
for stops in the dilepton + jets +$jt final state. Taken 
from [281. 
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FIG. 17: The mass of the CHAMP candidates as measured 
using the time-of-flight system. Taken from [29]. 
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IV. GAUGE MEDIATED SUSY BREAKING 
MODELS 



FIG. 18: The 95% confidence level cross section upper lim- 
its as a function of CHAMP mass in stop scenarios. Taken 
from [291. 



Gauge Mediated Supersymmetry Breaking (GMSB) 
models typically have the same mass relationships as in 
mSUGRA except for the G, the SUSY partner of the 
as- yet unobserved graviton. In mSUGRA models, the 
graviton is expected to have a mass at the GUT scale, 
but in GMSB models the G is expected to have a mass of 
a keV, making it the LSP. Thus, the decay X?— > 7G pro- 
vides a warm dark matter candidate consistent with as- 
tronomical observations and models of inflation [20]. It 
has the further advantage that it provides a more natural 
solution for FCNC problems than does mSUGRA. 

In typical models that take into account cosmological 
constraints [56] the Xi is favored to have a lifetime on the 
order of nanoseconds. This would make them stable on 
the timescale of the early universe, before they decayed 
into a G. Since limits on GMSB models exclude above 
myo < 100 GeV/c 2 [57], this predicts that squarks and 
gluinos are too massive to be produced at the Tevatron. 



So, as in mSUGRA, gaugino pair production dominates. 
Each gaugino typically produces a Xi> an d other high 
Et, light particles; the Xi decays via 7G. The G, like 
the Xi in mSUGRA models, leaves the detector and can 
give significant ftp. The lifetime and mass of the Xi dic- 
tate the different final states [58]. Specifically, if both 
Xi decay in the detector, a final state of 77 + pr + X is 
favored for a short lifetime (t^o < 1 ns). For intermediate 
lifetimes, 1 < t^o < 50 ns, the 7 + fr + X final state 
may be detectable if one Xi decays in the detector while 
the other leaves the detector without decaying or inter- 
acting. In this scenario, the time of arrival of the photon 
at the calorimeter can be the signal of a "long" lifetime 
Xi- For large lifetimes, both neutralinos can leave the 
detector and are indistinguishable from mSUGRA sce- 
narios. The different lifetime scenarios are considered 
separately. 
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A. Short-Lived 

For short lived Xi the dominant search channel is 77 
+ fir- CDF has recently developed a new tool to assess 
the significance of the fir measurement [34] . This allows 
for a straightforward separation of QCD backgrounds 
with no intrinsic fir from those with real fir such as 
EWK sources. Also a new photon timing device is used 
to remove non-collision sources of background [33]. A 
signature-based search in 2.6 fb _1 of data for 77 events 
with significant fir is shown in Figure 19 [30, 34]. To 
optimize the 77 + fir search for r%a <c 1 ns, the GMSB 
model analysis requires both significant fir and large Hr 
to identify the extra particles produced in the decay of 
the gauginos [30]. After optimization, the results are 
shown in the bottom of Figure 19. Since there is no evi- 
dence for new physics, cross section limits are set on spar- 
ticle production as a function of the Xi mass in Figure 20. 
Comparable results from D0 that assume r^o <C 1 ns can 
be found at [59] . Of particular note is that the null re- 
sults can also be interpreted as limits on production as 
a function of Xi lifetime as shown in the bottom of Fig- 
ure 20. 
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used to measure the time of arrival of photons in the de- 
tector [33], Figure 21 (Top) shows the technique for dis- 
tinguishing prompt photon production from "delayed" 
photons from long lifetime Xi that travel macroscopic 
distances then decay in the detector. The search for 
long-lived neutralinos requires a photon, large fir and 
a jet from the gaugino decays [31] in 570 pb _1 of data. 
Figure 21 (Bottom) shows the photon time of arrival cor- 
rected for the collision time and the time of flight for the 
sample. There is no evidence for new physics. Figure 22 
shows the 95% confidence level cross section upper limits 
in the lifetime vs. mass plane. Figure 23 shows the com- 
bined exclusion region from both analyses. Projections 
indicate that with 10 fb _1 the search sensitivity should 
be well into the cosmology favored region. 



FIG. 19: (Top) The fir significance distribution in inclusive 
77 events. (Bottom) The Ht distribution of 77 + fir events 
in the search for low lifetime Xi m GMSB events. Taken 
from [30]. 



CONCLUSIONS 



B. Long-Lived 

Cosmological constraints on GMSB favor keV 
G masses and nanosecond Xi lifetimes for xi masses 
above 100 GeV/c 2 . At CDF a photon timing system is 



The CDF experiment has performed a broad and deep 
set of cosmology inspired searches for SUSY. Unfortu- 
nately, it has found no evidence for new physics. As of 
mid-2009, the Tevatron is still delivering data at unprece- 
dented levels and the detectors are functioning beauti- 
fully. With data taking again started in 2009, the future 
is bright for discovery at the Tevatron. Until the LHC 
starts running, it is still leading the search for SUSY. 
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FIG. 21: (Top) A diagram of the how photons from long lived 
decays of Xi —* arrive at the detector with delayed times. 
(Bottom) The corrected timing distribution for the delayed 
photon analysis in the 7 + jet + ftrr search for long lifetime 
Xi GMSB models. Taken from [31]. 
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APPENDIX A: NOTES ON THE LHC STARTUP 

It is important to note that most of the analyses pre- 
sented in this paper will look like they were easy in hind- 
sight. At the time of this writing it's 2009 and Run II 
is eight years into into it's lifetime. These analyses are 
a lot more difficult than they look and take a lot longer 
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FIG. 22: Contours of constant 95% confidence level cross sec- 
tion upper limits on GMSB SUSY production from the 7 + jet 
+ pr final state search for long-lived Xi's. Taken from [31]. 
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FIG. 23: The combined GMSB exclusion regions in the life- 
time vs. mass plane of the x^- Taken from [30]. 



than they, perhaps, should. Run II at the Tevatron was 
effectively started in 2001, the first indirect SUSY search 
paper appeared in print in 2004 [43] with the first direct 
search paper following suit in 2005 [57]. We also note 
that the first direct search, in 77+^r, was the "sim- 
plest" search analysis with the photons dominating the 
calorimeter response, minimizing the fir resolution and 
the reconstruction pathologies, crucial for robust SUSY 
searches. Although not strictly true, both analyses were 
also largely repeats of Run I analyses [60, 61]. How 
quickly will the LHC and its detectors come up? CDF 
was a known detector and the Tevatron was a known 
accelerator, and there were five times fewer collaborat- 
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ing scientists. While there's only a sample of one and 
there will be great pressure to get results out quickly, a 
historical perspective says that since the CDF detector 
was known better the LHC results will come out more 
slowly relative to first collisions. The same time-ordering 
is to be expected since the Tevatron/CDF combination 



had working/known accelerator and well tested Monte 
Carlo simulations of the collision process and detector 
response. As for having five times as many people work- 
ing on CMS/ ATLAS than on CDF/D0, only time will 
tell if this makes things go quicker or more slowly, but 
compelling arguments can be made either way. 
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